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(g) Verfahren zur Bestimmung der Tern peraturvertei lung eines NOx-Katalysators 

@ In einem Verfahren zur Bestimmung der Tern peratur ei- 
nes Katalysators einer Brennkraftmaschine wird der Kata- 
lysator in n Scheiben, n > 1, in seiner axialen Richtung 
aufgeteilt und die Temperatur jeder Scheibe als Funktion 
der Temperatur des die Scheibe anstromenden Abgases 
bestlmmt, wobei die radiale Temperaturverteilung als 
konstant angenommen und eine adiabattsche Warme- 
ubertragung zwischen dem Abgas und der Katalysator- 
scheibe n berechnet wird. 
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Beschreibung 

rOOOn Die Erfindung betrifft ein Verfahrcn zur Bestinunung der Temperaturverteilung eines ^atalysators insbeson^ 

f»u~„\i^rTiMrielv!modusbeimagerbetreibbarenBrennkraftmaschinenbesUmm^ . ^ u 

[SS SfS^S-eSbTJe Brennto^^ beispielsweise dkekteinspritzende Brer^kraftm-chm^^^^ , 

[0002J "^"'^Semagerc^Lr „ . erzeugen wahrend des Magerbetriebs umweltschadigende Stickoxide 

N^'^SSS SeSe^t^.^^^^^^ der Brennkraf^aschine in ei.e. der Br^nkraftgaschine 

u* I u tT,^ MO^ Sneicherkatalvsator eespeichert. Aufgrund der endHchen Aufnahmekapazitat des NOx-Spei- 
TeSyS wiS d? SSSS'nS u'S gleichzeitiglr Reduktion zu Sdckstoff ^l^-^<^;; f^^J^^'^f:;^ 
?e JbeS to B^ennkrato freigesetzt Um zur Verminderung des Kraftstoffverbrauchs und der Emissionen den 

Fettbemeb . • / mogUch steuem zu kQnnen. sind moglichst genaue Kenntnisse der 

ISeTmAbtiS^ 

KasSSe gender auf die Randbedingungen und Grenzen des Abgasreinigungssystems abzusUmmen. wobei msbe. 
sondere zu nennen sind: 

- Anreicherung eines atomotors im oberen TeiUast- und VoUastbereich zur Begrenzung der Katalysatortempera- 

"iTu^^n'^rab'gt'Soder katalysatortexnperatursteigemden MaBnahmen nach AuskUhlung und nachfol- 
eendem Erwarmen eines hinreichend groBen Bereichs des Katalysators und k^j a„ Vnt- 

! sSSlen einer hinreichenden &warmung eines ausreichend groBen Bereichs des Katalysators bei der Ent- 

schwefelung. 

r00041 Beim gegenwartigen Stand der Technik kann die Katalysatortemperatur beispielsweise durch ^J^^en der Ab- 
L WerXsowohl vor L auch hinter dem Katalysator ermittelt werden. Sie kann auch durch lessen der Abg^^^^^ 

B^enntaftmaschine richten sich danh besser nach dem realen Verhalten des Abgasreimgungssystems, da die Emgnffs- 
schwellen naher an die kritischen GroBen des Katalysators herangefiihrt werden konnen. R^^^hnnno 

P^:^S^uliS.i£^^i^^n <1M uber den MogUchkeiten gegenwartiger und in naher Zukunft zu erwartender 

SrlSung liegt daher die Aufgabe zugrunde. ein Verfahren zur Berechnung der Temperaturverteilung zu 
i?J!ffln das bd we"entbch geringerem Umfang und unter Inkaufnahme einer gewissen Ungenauigke.t Aussagen uber 
7MctT:m^Z^elSrn Katalysa.of em.ogUcht. Femer wird ein Verfahren zur Steuerung «nes Speicherka- 

. iSri^Au^a^vf^T^^^^^^^^^^ 

KysteiS™!e radiale TempeTaLverteilung als konstant angenommen, der Katalysator wud in n hintereinander ange- 
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ordnete Segmente unterteilt und die Temperatur jedes Segment bestimmt. Dabei ist n eine naturliche Zahl groBer Eins. 
Die jedem n zugeordneten Segmente konnen dabei unterschiedlich lang sein, um eine bessere Anpassung zu erreichen. 
Durch Unterteilung des Katalysators in n hintereinander angeordnete Segmente wird die Zahl der zu berechnenden Ka- 
talysatorelemente auf ein praxisgerechtes MaB begrenzt. 

[0010] Vorzugsweise wird der Katalysator in n Scheiben zu je 1/n seiner Gesamtlange unterteilt Die Zahl n kann dabei 5 
zwischen 1 und unendlich liegen; unter realen Berechnungsbedingungen wLrd n zwischen 2 und 50, vorteilhaft zwischen 
3 und 12 und insbesondere zwischen 5 und 7 betragen. 

[0011] Ein Warmeiibergang zwischen den Scheiben findet im wesentlichen nur durch das durchstromende Gas statt. 
Die axiale Warmeleitung wird daher vorzugsweise z,u Null gesetzt; sie kann jedoch fiir besondere Anwendungsfalle in 
die nachfolgend aufgeHsteten Berechnungsalgorithmen mit aufgenommen werden. lO 
[0012] Zunachst wird der adiabate Warmeiibergang zwischen dem Abgas und der Katalysatorscheibe und damit die 

durch thermische Energieiibertragung erzielte Temperaturanderung des Abgases und der Katalysatorscheibe berechnet. 
AnschUeBend wird abhangig vom Abgasmassenstrom und von der Temperatur jeder einzelnen Katalysatorscheibe eine 
ortlich variable Umsetzung der anstromenden Schadstoffe angenommen. Die Warmeverluste iiber den Katalysator wer- 
den naherungsweise in Abhangigkeit von der Temperatur der Katalysatorscheibe beriicksichtigt. is 
[0013] Zur Berechnung des thermischen Warmeiibergangs zwischen Abgas und Katalysator an jeder einzelnen Kata- 
lysatorscheibe, des Energieubergangs durch Schadstoffumsetzung und der Warmeverluste im Vorrohr konnen besonders 
vorteilhaft die im Anhang 1 und 2 beschriebenen Algorithmen genutzt werden. 

[0014] Die Abgastemperatur im Vorlauf des Katalysators wird durch einen Temperaturfiihler vor dem Katalysator oder 
durch einen Temperaturfiihler hinter einem eventuell vorgeschalteten Vorkatalysator ermittelt. Im letzten Fall muB der 20 
Temperaturverlust iiber das Abgasrohr zwischen Vor- und Hauptkatalysator modelliert werden. 

[0015] Zur Berechnung der Schadstoffumsetzung und der resultierenden Temperaturerhohung wird zunachst abhangig 
von der Raumgeschwindigkeit ein Profil erstellt, das den Umsetzungsanteil der Restschadstoffe an den einzelnen Kata- 
lysatorscheiben angibt. Wahrend bei niedrigen Raumgeschwindigkeiten die Reaktion iiberwiegend in der oder den vor- 
dersten Katalysatorscheiben stattfindet, verlagert sich mit zunehmender Raumgeschwindigkeit die Reaktion in die weiter 25 
hinten liegenden Katalysatorscheiben. Bei sehr hohen Raumgeschwindigkeiten kann sogar der Katalysator ganz oder 
teilweise tiberlaufen werden. Mit der Temperatur der aktuell berechneten Katalysatorscheibe wird aus einem Kennfeld 
die fiir diese Berechnungsschleife zutreffende Konvertierungsrate separat fur jeden einzelnen Schadstoff ermittelt. 
[0016] Die HC-, CO-, NOx-, 02- und H2- Rohemissionen konnen iiblicherweise nicht gemessen werden. Es ist jedoch 
eine naherungsweise Bestinunung aus einem Kennfeld mogUch, das aus den stationaren, betriebswarmen Emissionen 30 
mit additiver oder multiplikativer Korrektur fiir dynamische Vorgange und Kaltstarteffekte die realen Emissionen errech- 
net. Das Kennfeld wird in bekannter Weise auf einem Motorenpriifstand bestimmt. 

[0017] Da modeme Abgasreinigungssysteme zunehmend aus einem motomahen kleinvolumigen Vorkatalysator und 
einem nachgesichalteten Hauptkatalysator bestehen, wird der Hauptkatalysator iiberwiegend nicht mehr mit der gesamten 
Rohemissionsmenge beaufschlagt, sondem lediglich mit einer Teilmraige. Abhanjgig von der Raumgeschwindigkeit des 35 
Vorkatalysators wird ein Rohemissionsanteil ermittelt, der nicht vom Vorkatalysator umgesetzt wird und somit den 
Hauptkatalysator beaufschlagt. 

[0018] Das temperaturabhangige Konvertierungsverhalten des \brkatalysators nach einem Kaltstart muB dabei nicht 
beriicksichtigt werden, da davon auszugehen ist, daB der Vorkatalysator aufgrund seiner Anordnung, Beschichtung und 
Masse schnelier anspringt als der Hauptkatalysator und daB bei Einsetzen katalytischer Reaktionen auf dem Hauptkata- 40 
lysator der Vorkatalysator seine Arbeitstemperatur und damit seine maximale temperaturabhangige Konvertierungsrate 
bereits erreicht hat. 

[0019] Mit der Temperatur und der Raumgeschwindigkeit wird - wie oben beschrieben - die Konvertierungsrate der 
Katalysatorscheibe berechnet. Die anstromenden berechneten Rest-S chads toffe werden auf der Katalysatorscheibe ganz 
oder teilweise unter Beachtung des Restsauerstoffahgebotes umgesetzt. Die Energiebilanz an der Katalysatorscheibe ist 45 
anschUeBend nochmals unter Einbeziehung der umgesetzten chemischen Energie durchzurechnen, um die abstromende 
adiabate Abgastemperatur zu ermitteln. Die Abgas- und Katalysatorscheibentemperatur werden abschlieBend noch um 
den Warmeverlust korrigiert. 

[0020] Eine bevorzugte Ausfiihrungsform der Erfindung wird nachfolgend anhand der einzigen Figur erlautert, 
[0021] Die Figur zeigt schematisch den Ablaufplan der Berechnung des Temperaturprofils eines Katalysators fiir einen 50 
Zeitschritt wobei als Berechnungsfrequenz ist 0,1 bis 10 Hz, vorteilhaft um 0,5 bis 2 Hz und insbesondere sind 1 Hz 
empfehlenswert ist. 

[0022] Zur Berechnung der Temperaturverteilung des Katalysators werden verschiedene EingangsgroBen benotigt, die 
in zwei Nebenrechnungen ermittelt werden, wobei die EingangsgroBen zur Berechnung einer Scheibe i von den Ergeb- 
nissen der Berechnung der Temperatur der Scheibe i-1 abhangen konnen. 55 
[0023] Aus der GroBe Rohemissionen 1, die einem auf einem Priifstand ermittelten Kennfeld entnommen werden, so- 
wie einem Mot'ordynamik-Korrekturfaktor 2 und einem Kiihlmittel-Korrekturtaktor 3 werden in 4 die Emissionen vor 
dem Vorkatalysator berechnet. Die genannten Korrekturfaktoren 2 und 3 werden auf bekannte Weise bestimmt. In 5 wird 
aus den in 4 errechneten Emissionen die katalytische Umsetzung der Emissionen im Vorkatalysator berechnet, wobei zu 
dieser Berechnung eine weitere EingangsgroBe notwendig ist, und zwar die im Schritt 7 bestimmte Konvertierungsrate 60 
des Vorkatalysators fiir jeden Schadstoff, die abhangig ist von der Raumgeschwindigkeit des Abgases. Die Raumge- 
schwindigkeit des Abgases wird vom Block 6 bereitgestellt. Aus der im Block 5 erfolgten Berechnung der katalytischen 
Umsetzung der in 4 korrigierten Rohemissionen 1 ergibt sich die vor dem Katalysator anstehenden Emissionen 8, ge- 
nauer gesagt die Emissionen 8, die vor der ersten Katalysatorscheibe anstehen. Bei diesen Emissionen handelt es sich um 
die HC-, CO-, NOx-, H2 und 02-Massenstrome. 65 
[0024] Mittels einer Messung der Luftmasse in Schritt 9 und der Messung der Temperatur in Schritt 10 wird in Schritt 
11 der Massenstrom und die Temperatur des Abgases hinter dem Vorkatalysator bestimmt. Mittels beispielsweise des im 
Anhang 2. dargestellten Algorithmus wird im Schritt 12 eine Berechnung der Einfliisse des Vorrohres vorgenommen und 
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SS^^r^^runfer^^^^^^ von der n-ten Scheibe abstr.n.enden Abgases 

t Abg.n+1 bestimmt 

Anhang 

[0027] Fiir die hier anschUeBenden Betrachlungen werden die folgenden Definitionen getroffen: 
AKat Katalysatoroberflache 
AKRohr Oberflache Katalysatorcanning (Hiillrohr) 
AroHt Rohroberflache 

Cp^bg spezifische Warmekapazitat des Abgases , 
40 Cp Kat spezifische Warmekapazitat des Katalysators, 

Cp Rohr spezifische Warmekapazitat des Rohres 

Hu,schad5toff Heizwert der Schadstoffe 

mi^ Katalysatormasse, 

niRohrRohnn^se 
45 ihschadstoffSchadstoffmassenstrom ^ 

riiAbg Abgasmassenstrom 

Rk Konvertierungsrate 

tAbg Abgastemperatur 

tKat Katalysatortemperatur 
50 t'Kat Katalysatortemperatur am Eintritt in die Katalysatorscheibe 

tu Umgebungstemperatur 

tRohi Rohrtemperatur 

a Warmeiibergangskoeffizient 

ttAK Warmeiibergangskoeffizient Abgas-Katalysator 
55 ttAR Warmeubergangskoeffizient Abgas-Rohr 

ttKRU Warmeubergangskoeffizient Katalysatorcannmg-Umgebung 

At Zeitintervall 

T Zeit 

[0028] Weiter haben die verwendeten Indizes folgende Bedeutung: 
60 ' Eintritt in die Katalysatorscheibe 
i ortlicher Laufindex 
T zeitlicher Laufindex 
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Berechnungsschleife fiir eine Katalysatorscheibe 

Berechnung der thermisch adiabatischen Energieiibertragung zwischen dem Abgas und der betrachteten Katalysator- 
scheibe zur Ermittlung der Konvertierungsrate der Katalysatorscheibe 

[0029] AUgemein laBt sich die Erwarmung eines Katalysators durch die folgende Dififerentialgleichung (1) beschrei- 
ben: 

'"^Ca. C^JCa,'^=CC^ A^^(t'^^-t'^{T)). (1) 

[0030] Dabei erfullt die obige Differentialgleichung (1) folgende Randbedingung: 

tWi:=0) = Tu (2) . • - 

[0031] Aus der formalen Losung der obigen Differentialgleichung unter Beriicksichtigung der Randbedingung laBt 
sich die Temperatur des Katalysators in ortlichen (Laufindex i) und zeitlichen Schleifen (Laufindex T) enechnen. Als Er- 
warmung ergibt sich am Eintritt in die Katalysatorscheibe i: . 


'1 


[0032] Bei einer zeitlichen Anderung der Abgastemperatur wird zur Berechnung der Laufindex X neu gestartet und als 
Randbedingung die bisherige Katalysatortentiperatur tKat^x-i eingesetzt. Es ist jedoch auch mdglich, daB mit dem Zeitin- 
tervall Ax und der Randbedingung fur die vorherige Katalysatortemperatur gerechnet wird. Der Warmeiibeigangskoeffi- 
zient a wird mit Hilfe der bekannten Nusseltzahl fiir den vorliegenden Stromungszustand naherungsweise ennittelt. 

Berechnung der Energieiibertragung durch die Schadstofifiimsetzung 

[0033] Mit Hilfe einer vorlaufigen Katalysatortemperatur kann eine Konvertierungsrate Rk der aktuellen Katalysator- 
scheibe aus einem Kennfeld mit den Eingangsdaten t'KaU und der ortlichen Raumgeschwindigkeit ermittelt werden. Bei 
der Berechnung der Energieiibertragung durch die Schadstoffumsetzimg im Katalysator wird sowohl der voirhandene 
Restsauerstoffgehalt als auch die Konvertierung im Vorkatalysator beriicksichtigt. Aus der Energiebilanz laBt sich in be- . 
kannter Weise die Temperatur des Abgases berechnen: 

f __ '^Abg '^p^g '^AbgJ '^Rk '^Scha dstoff ' ^ajSchadstoff 

^Abg4- -r (4) 

^Abg'^p^g 

[0034] Aus dieser Abgastemperatur, die die chemischen Umsatze beriicksichtigt, laBt sich die reale Temperatur der be- 
trachteten Katalysatorscheibe wie folgt ermitteln: ' 


^Kmj (^) = ^AtgJ t (^ICat,i,r^l " Ugj ) ' CXp 


(S) 


[0035] Durch die Warmeabgabe des Abgases an die Katalysatorscheibe kiihlt sich das Abgas ab und aus der Energie- 
bilanz 

^Abg • ^p^bg ' ^AbgJ = ^Abg ' <^p^Abg ' ^Abg^M ^AK " ^Kat ' i^Abg^i ^ ^ Kai,i ) (6) 

ergibt sich die Abgastemperatur an der nachsten Katalysatorscheibe: . 

- ^Abg '^p^Abg '^Abg,i ''^AK'^Kat 'ihbgj ^^Kaij) _ 

Ubrg^^l (7) 

^Abg '^p,Abg 

Beriicksichtigung von Verlusten 


[0036] An der Katalysatoroberflache treten Warmeverlust, beispielsweise durch Abstrahlung oder Konvektion, auf, so 
daB sich ein radialer Temperaturgradient ausbildet. Diese Warmeverlust werden iiberschlagsmafiig bestimmt. Wederum 
iiber eine Energiebilanz 
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wild dann die sich ergebende Temperatur der Katalysatorscheibe ermittelt: 

Berechnung des Vorrohrs 

[0037] . Die Erwarmung des Vorrohrs, das als isotherm angenommen wird. laBt sich wie folgt berechnen: 
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25 


(10) 


folgen, wobei in dieser Gleichung die OroBen des ^^^^^f " ^^bei diese \ferluste durch Konvektion 

vSnsverluL. erjbt sich schUeBUch eine Teinperamr des Vorrohrs von: 


'1 


■RU '^Rahr (11) 


BEZUGSZEICHENLISTE 

1 Rohemissionen aus Kennfeld 

2 Korrekturf aktor Motordynamik 

35 3 Korrekturf aktor Kuhlmitteltemperatur 

4 Emissionen vor dem Vorkatalysator 

5 Katalytische Umsetzung im Vorkatalysator 

6 Raumgeschwindigkeit des Abgases 

7 Konvertierungsrate des Vorkatalysators bei Arbeitstemperatur 
40 8 Emissionen vor der ersten Katalysatorscheibe 

9 Luftmassenmessung 

10 Temperaturmessung Abgas nach Vorkatalysator 

11 Abgas nach Vorkatalysator 

12 Berechnung Vorrohr 

45 13 Abgas vor erster Katalysatorscheibe 

14 vorlaufigeTemperatur der ersten Katalysatorscheibe 

15 Temperaturerhohung durch Exothemiie der ersten Katalysatorscheibe 

16 Konvertierungsraten der ersten Katalysatorscheibe 

17 Abgaskuhlung 

18 Emissionen vor der zweiten Katalysatorscheibe 

19 vorlaufige Temperatur der zweiten Katalysatorscheibe 

20 Temperaturerh5hung durch Exothermic der zweiten Katalysatorscheibe 

21 Konvertierungsraten der zweiten Katalysatorscheibe 

22 Abgaskuhlung 

23 Emissionen vor der n-ten Katalysatorscheibe 

24 vorlaufige Temperatur der n-ten Katalysatorscheibe 

25 Temperaturerhohung durch Exothermie der n-ten Katalysatorscheibe 

26 Konvertierungsraten der n-ten Katalysatorscheibe 

27 Abgaskuhlung 

Patentanspriiche 
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Katalysators aufweist. 

3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB die Scheiben unterschiedliche Langen aufweisen 

4. Verfahren nach einem der vorangegangenen Anspriiche, dadurch gekennzeichnet, daB fiir n gilt: 2 < n < 50. 

5. Verfahren nach Anspruch 4, dadurch gekennzeichnet, daB fur n gilt: 3 < n < 12. 

6. Verfahren nach Anspruch 4, dadurch gekennzeichnet, daB fiir n gilt: 5 < n < 7. 5 

7. Verfahren nach einem der vorangegangenen Anspriiche, dadurch gekennzeichnet, daB die axiale Warmeleitung 
gleich Null gesetzt wird. 

8. Verfahren nach einem der vorangegangenen Anspriiche, dadurch gekennzeichnet, daB die Temperaturanderung 
einer Katalysatorscheibe durch chemische Umsetzung der eine Katalysatorscheibe n anstromenden SchadstofFe als 
ortlich variable Funktion der zuletzt berechneten Katalysatorscheibentemperatur und des Abgasmassenstroms aus lo 
der Energiebilanz berechnet wird. 

9. Verfahren nach Anspruch 8, dadurch gekennzeichnet, daB die Temperaturanderung einer Katalysatorscheibe n 
durch Warmeverluste iiber den Katalysator als Funktion der zuletzt berechneten Temperatur der Katalysatorscheibe 
n bestinunt wird. . 

10. Verfahren nach Anspruch 9, dadurch gekennzeichnet, daB der Umsetzungsanteil der Schadstoffe fiir jede Kata- 15 
lysatorscheibe als Funktion der Raumgeschwindigkeit bestimmt wird. 

1 1 . Verfahren nach Anspruch 10, dadurch gekennzeichnet, daB fiir jeden SchadstoflF separat die Umsetzungsrate be- 
stimmt wird. 

12. Verfahren nach Anspruch 11, dadurch gekennzeichnet, daB die Umsetzungsraten der Schadstoffe einem Kenn- 
feld entnommen werden. 20 

13. Verfahren nach einem der vorangegangenen Anspriiche, dadurch gekennzeichnet, daB die Berechnungsfre- 
quenz zur Berechnung eines Temperaturprofils zwischen 0,1 bis 10 Hz liegt. 

14. Verfahren nach einem der vorangegangenen Anspriiche, dadurch gekennzeichnet, daB die Abgastemperatur im 
Vorlauf des Katalysators durch einen Temperaturfuhler ermittelt wird, wobei der Temperaturverlust des Abgases 
iiber das Abgasrohr zwischen dem Ort des Fiihlers und dem Eintritt des Abgases in den Katalysator bei der Berech- 25 
nung beriicksichtigt wird. 

15. Verwendung des Verfahren nach einem der vorangegangenen Anspriiche zur Steuerung des Betriebsmodus ei- 
ner magerbetreibbaren Brennkraftmaschine in eihem Motorsteuergerat zur Bestimmung des Betriebsmodus der 
Brennkraftmaschine. 

30 

Hierzu 1 Seite(n) Zeichnungen 
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